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We investigate nonequilibrium underdamped Langevin dynamics of Brownian particles that in¬ 
teract through a harmonic potential with coupling constant K and are in thermal contact with two 
heat baths at different temperatures. The system is characterized by a net heat flow and an entropy 
production in the steady state. We compare the entropy production of the harmonic system with 
that of Brownian particles linked with a rigid rod. The harmonic system may be expected to reduce 
to the rigid rod system in the infinite K limit. However, we find that the harmonic system in the 
A —>■ 00 limit produces more entropy than the rigid rod system. The harmonic system has the center 
of mass coordinate as a slow variable and the relative coordinate as a fast variable. By identifying 
the contributions of the degrees of freedom to the total entropy production, we show that the hidden 
entropy production by the fast variable is responsible for the extra entropy production. We discuss 
the K dependence of each contribution. 

PACS numbers: 05.70.Ln, 05.40.-a, 02.50.-r, 05.10.Gg 


I. INTRODUCTION 

Recent developments in stochastic thermodynamics 
boost active researches on microscopic systems in ther¬ 
mal environments. Stochastic thermodynamics allows 
one to define thermodynamic quantities such as heat and 
work at the microscopic trajectory level, with which ther¬ 
modynamic properties are readily studied 0,i. 

When one studies a thermodynamic system, a coarse¬ 
grained description may be necessary [3l4l2l|. Suppose 
that some transition rates among microscopic configura¬ 
tions are much faster than the others. Such systems can 
be coarse-grained in terms of so-called mesostates which 
refer to a set of microscopic configurations connected via 
fast processes The coarse-graining is also necessary 

if some degrees of freedom are inaccessible due to e.g., 
a resolution limit of measurement devices I, [13. The 
coarse-graining naturally raises questions on the formu¬ 
lation of effective thermodynamics and on the extent to 
which the thermodynamic laws are valid. 

In this paper, we investigate the entropy production 
by a harmonic chain of Brownian particles between two 
heat baths at different temperatures with the focus on 
the the role of fast and slow degrees of freedom. When 
two heat baths are connected with a medium, heat flows 
from a higher temperature bath to a lower temperature 
one and the total entropy increases [iM3- A toy model 
for this phenomenon in one dimension was studied thor¬ 
oughly in Refs. EH with the equation of motion for 
the displacement X: 

dX ^ , dX ^ , 

where ji (i = 1,2) is the damping coefficient of a 
heat bath i at temperature Ti and ^i{t) is the thermal 
Langevin noise satisfying 

im) = 0 , = 2^,TAjS{t-1') . ( 2 ) 


Ti T\ K ^2 

(a) (b) 

FIG. 1. Comparison of a rigid rod system in (a) and a har¬ 
monic system in (b). 
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The Boltzmann constant fcs is set to unity. The model 
describes a single Brownian particle of mass M coupled 
to two heat baths simultaneously. At the same time, as 
illustrated in Fig. [Ua), it also describes two Brownian 
particles of total mass M tied with a ri gid rod. This 
rigid rod system was also studied in Refs. flO - H^ for the 
heat conduction in the Feynman ratchet |19| . 

We notice that the model with a single degree of free¬ 
dom assumes an ideal medium for heat conduction, which 
could be realized with a many-body system in the strong 
coupling limit. A many-body system has multiple degrees 
of freedom such as the center of mass coordinate and the 
relative coordinates. In the infinite coupling limit, the 
relative coordinates relaxing fast can be coarse-grained 
out to leave the center of mass coordinate as a relevant 
slow variable. It is interesting to investigate the simi¬ 
larity and the dissimilarity between the system given by 
(P) and a system with multiple degrees of freedom in the 
strong coupling limit. 

For the purpose, we consider two Brownian particles 
exchanging a harmonic interaction as shown in Fig.pb). 
Heat conduction and the entropy production in harmonic 
systems have long been studied [ij, [20l - [^ . In contrast 
to the previous works, our study focuses on the role of the 
slow and fast variable. Harmonic systems are analytically 
tractable [2^-|3l|. We derive the formula for the entropy 
productions by the slow and fast variables as a function 
of the coupling constant, which will shed a light on the 
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effect of the coarse graining. 

This paper is organized as follows: The model sys¬ 
tem consisting of two Brownian particles is introduced 
in Sec. m The heat flow and the entropy production in 
this model are calculated in Sec. mil The slow and fast 
variable are introduced and their contributions to the to¬ 
tal entropy production are investigated in Sec. IIVI We 
compare the harmonic system and the rigid rod system, 
and draw the conclusion for the hidden entropy produc¬ 
tion by the fast variable. We conclude the paper with 
summary and discussions in Sec. El 


II. TWO PARTICLE SYSTEM 


Two particles of mass m in one dimension interact 
through a harmonic potential 

V{xi,X2) = ( 3 ) 

where Xi {i = 1,2) is the displacement of a particle i and 
K is the coupling constant. A particle i is in contact with 
the ith heat bath which is characterized with a tempera¬ 
ture Ti and a damping coefficient 7 ^ (see Fig.jljb)). The 
underdamped Langevin equation is given by 


Xl = Vl, 

X2=V2, 

mvi = -71U1 - K{xi- X2) 

mv2 = -J 2 V 2 - K{x 2 - Xl) - 1 - 6 ( 0 ) 

where { 6 ( 0 } represent the thermal noises satisfying ([ 2 |). 

Suppose that the system evolves from a configura¬ 
tion (xi,vi;x 2 ,V 2 ) at time t to {xi + dxi,vi dvi ; X 2 + 
dx 2 ,V 2 + dv 2 ) in an infinitesimal time step dt. Accord¬ 
ing to stochastic thermodynamics [I|, the total entropy 
production is given by 

dStot = dSs + dSb, (5) 


where dSs denotes the Shannon entropy change of the 
Brownian particles and dSb is the entropy change of the 
thermal heat baths. The latter is given by 


dSb 


dQi dQ2 


( 6 ) 


where 


dQi = ViO {-^iVidt + dWi) (7) 


is the heat absorbed by the Brownian particle from the 
bath i, and dWi = dt'^i(t')dt' satisfies (dWi) = 0 
and {{dWi)'^) = 2^iTidt. The notation o means that the 
product is in the Stratonovich sense [3^ . 

The heat can be expressed in different ways. By re¬ 
placing {—^iVidt + dWi) with {mdvi+K{xi—X 2 )dt), we 


obtain 


/I 2 1 2 1 

dQi =d ( -^mvi + -^Kxi - -KxiX2 

+ —{X 1 V 2 - X2Vi)dt. 


Similarly, we also obtain 
'1 


dQ2 =d 


1 


-b ^Kx 2 - ^KxiX2 

K, 


- —{xiV 2 - ViX 2 )dt. 


( 8 ) 


( 9 ) 


They add up to dE =dQi +dQ 2 , where dE is the change 
in the the total internal energy E = ^m{v‘f -b u|) -b 
V{xi,X 2 )- The total derivatives in (jS]) and ([S]) do not con¬ 
tribute to the steady state average. Hence, the steady- 
state average of the heat flux, denoted by qi = {dQi/dt)s 
is given by 

K 

q = qi = -q 2 =-j-{xiV 2 - X 2 Vi)s , (10) 

where (•)s means the steady state ensemble average. 

The heat flux acquires an interesting interpretation by 
regarding Xi and Vi as the coordinates of the position vec¬ 
tor X = {xi,X2,0) and the velocity vector v = (ui,t;2,0) 
of an imaginary particle moving on a xy plane in a three 
dimensional space. Then the heat flux in (HOD becomes 
proportional to the z component of the angular momen¬ 
tum L = mx X D of the imaginary particle: 

g=^(A3). . ( 11 ) 

A nonzero heat flow implies a rotational motion. Such 
a rotational motion induced by a nonequilibrium driv¬ 
ing was investigated in the context of a Brownian gyra- 
tor [13 ■ 

We find another useful expression for the heat flow. 
The Stratonovich product vi o dWi in 0 is defined as 
i[ui-b(ui-bdui)]dWi dll , where dui = [—'yiVidt—K{xi — 
X 2 )dt+dWi]/m from the Langevin equation. The ensem¬ 
ble average yields (viodWi) = Thus, the average 

heat flux in the steady state can be written as 


9 


hi (h 

m \ 2 



( 12 ) 


The heat flow is driven by the imbalance between the 
temperature of the heat bath and the effective tempera¬ 
ture of the particle measured by the kinetic energy. Note 
that the expression (113 is valid not only in the steady 
state but also in a transient state. 

In the steady state, the Shannon entropy of the par¬ 
ticles does not change {{dSs)s = 0). Consequently, the 
total entropy production rate in the steady state is given 

by 





(13) 


The heat flux and the entropy production rate is propor¬ 
tional to each other. 
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III. HEAT FLOW AND ENTROPY 
PRODUCTION 


For the heat flux, we need to evaluate the second mo¬ 
ments of the velocities, My = {vi{t)vj{t))s. The model 
with the linear Langevin equation in Q belongs to the 
class of the multivariate Ornstein-Uhlenbeck process . 
The second moments of such a system satisfy a set of 
linear equations , from which My can be easily eval¬ 
uated. Instead of using the method, we evaluate the mo¬ 
ments explicitly for self-containedness. 

The Langevin equation is solved in the Fourier 
space. We introduce column vectors x{uj) = 
(ii(w),i2(a;))'’’ and |(w) = (^i(a;), |2(a;))'’’ with Xi{uj) = 
jToo and |i(a;) = The super¬ 

script denotes the transpose. The Fourier components 
of the thermal noise obey 

= 0 ) . ( 14 ) 


The Langevin equation becomes F(u;) 5 (w) = ^(w), where 
F(w) is a 2 X 2 matrix given by 

p/ , _ /—— ijiUJ + K —K \ 

—K —mu}"^ — 1720; + KJ ' 


The steady state solution is given by x { uj ) = F(w)“^^(a;). 
The homogeneous solution is ignored because it vanishes 
in the steady state. For convenience, we write the solu¬ 
tion as x(uj) = where 

D(uj) = (—i7iw-FA")(—J72W-FAr)—( 15 ) 

is the determinant of F(w) and C(u;) = C(a;)"'" is the co¬ 
factor matrix given by 

p/ s _ f+ K K \ 

K —mux^ — i'jiuj + KJ ' 

( 16 ) 

By using Vi{u!) = —iu}Xi{uj) and the noise correlations 
in m , we can write 


for the heat flux and 

7i72i^ (Ti - T2)^ 

ill + l 2 )iKrn + ^1^2) T1T2 


( 20 ) 


for the entropy production rate. 

The heat flow is nonzero when Ti T2 and K is 
nonzero. The system absorbs a heat from a higher tem¬ 
perature heat bath and dissipates to a lower temperature 
one, which increases the entropy of the whole system. 


IV. HIDDEN ENTROPY PRODUCTION 


The heat flux and the entropy production rate increase 
monotonically with the coupling constant K. In the K —> 
00 limit, the heat flow rate converges to 


qoo = lim q 

K—^oo 


7172 

to( 7 i + 72) ^ ^ 


T2), 


( 21 ) 


and the entropy production rate converges to 


7i72 (^1 - '^ 2 )^ 

^(71+72) T1T2 


( 22 ) 


In this limit, the relative displacement of the two parti¬ 
cles (1x2 —Xi I ^ vanishes. Hence, one may regard 

the limiting case as the rigid rod system, that is, a single 
Brownian particle of mass M = 2 m (see Fig. [ij. It is 
worthy to compare the heat flux in both cases. The heat 
flux gi-igid of the rigid rod system is given by Ref. [H, [l^ 


Q'rigid — 


7i72 


M(7i -F 72) 


{Ti - T 2 ) = 


7i72 


2m (71 -F 72) 


(Ti - T 2 ). 


(23) 

Surprisingly, the heat flux frigid is the half of the limiting 
values (7oo because of the factor 2 in the denominator. 

In order to understand the origin of the discrepancy, we 
rewrite the Langevin equation (jT]) in terms of the center 
of mass coordinate X = {xi + X2)/2 and the relative 
coordinate x = xi — X 2 - The equations of motion for X 
and x are given by 


My = 


duj liTiCgiu})Cji{-u;)] 

TT D{uj)D{—uj) 


( 17 ) 


This integral is evaluated using the contour integral 
method, as detailed in Appendix |Al The results are 


Mn = -- - ^2) 

m (71 + 72)(Aim-F 7172)’ 


M 22 = 


T. 


71 A(ri - T2) 


m (71-F 72)(Arm-F 7172) ’ 
Mi 2 = M 21 = 0 . 

Inserting these into (TT^ . we finally obtain 
7 i72Ar(Ti - T 2 ) 


q = 


(71 +72)(Arm-F7i72) 


( 18 ) 


( 19 ) 


A = U, 

X = n, 

Mr = -(7i + 72)r-hl.^„ + d((), 

where M = 2 m is the total mass, /r = m /2 is the reduced 
mass. Cl = Cl+6 obeying (Ci(t)) = 0 and (Ci( 0 Ci(^')) = 
2(7iri -F 72T2)^(t - t'), and C2 = (Ci - C2)/2 obeying 
(C2(0) = 0 and (C2(i)C2(i')) = ^71^1 +-/2T2)S(t - t'). 
The center of mass and the the relative coordinates 
are coupled through the noise correlation {C,i{i)C,2{i')) = 
(71F1 — "i 2 T 2 ) 5 {t — t') and the friction forces. 

We can grasp an important feature of the system by 
considering the case with 71 =72 =7. The center of 
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mass performs a Brownian motion with the effective tem¬ 
perature Tgg = (Ti -I- T 2 )I 2 , and the relative coordinate 
performs the Ornstein-Uhlenbeck process under the 
harmonic potential with the same effective tem¬ 

perature Teff. The center of mass coordinate diffuse with 
the time scale tx ~ m/7. The relative coordinate os¬ 
cillates with the time scale Tx yJrnjK with the am¬ 
plitude I a; I ~ yjTes/K. Therefore, in the K ^ 00 
limit, X becomes a fast variable {tx tx) and becomes 
frozen (|x| —^ 0 ) in amplitude, while X becomes a slow 
variable being governed by the same equation of motion 
in m as the rigid rod system. 

Although the fast variable x becomes frozen in the 
K ^ 00 limit, it may leave a signature in thermodynamic 
quantities. We scrutinize this possibility by decompos¬ 
ing the heat flux, or equivalently the entropy production 
rate, into the sum of contributions of the fast and the 
slow variables. We start with the expression (|S]) for gen¬ 
eral 7i and 72. After inserting vi = V + v /2 into ([8]), we 
obtain 

^Qi ={—"fiV^dt + Vo dWi) + {—"fivVdt) 

+ 


Each term corresponds to the heat flow via the slow vari¬ 
able, the interplay between the slow and fast variables, 
and the fast variable, respectively. Thus, the average 
heat flow rate in the steady state is decomposed as 

q = qs + QF + Qm , ( 26 ) 

where 

« = ^ (t - 


is the contribution from the slow variable. 





is from the fast variable, and 


( 28 ) 


qm. = -li{vV)s ( 29 ) 


is the mixed contribution. 

Each contribution is easily evaluated using the second 
moments obtained in m- For example, {V^)s = 
(Mil + 2 Mi 2 -I- M22)/4. The results are 


qs = 


7172 2 Km -b 7i -b 7172 


Qm — 


4 m (71-b 72)(A'm-b 7172) 
7 i 7 i 72 


2 m {Km + 7172) 


qp = qs 


(Ti - T2) ( 30 ) 

(Ti - T2) ( 31 ) 

( 32 ) 



FIG. 2 . Heat flow rates as a function of K. Parameter values 
are 71 = 2 , 72 = 1 , m = 1 , and {Ti — T2) — 1 . 


mixed contribution q^n has the opposite sign to qs and 
qp, and vanishes in the infinite coupling limit. They are 
plotted as a function of K in Fig. [ 2 ] 

Measurements are done within the resolution limit of 
probing devices. Suppose that one experiments on the 
harmonic system with large K with a device whose time 
resolution T^es ^ Fx- Then, the harmonic system would 
be observed as the rigid rod system. Measurements 
would yield the apparent entropy production from the 
slow variable only, which is a part of the total entropy 
production. 

The opposite case is also interesting. When the cou¬ 
pling constant vanishes {K —>■ 0 ), the two particles are 
in their own thermal equilibrium and the heat flux and 
the entropy production vanish. However, the apparent 
heat flux qs associated with the center of mass coordi¬ 
nate still remains finite even at iF = 0 . Suppose that the 
two noninteracting particles are confined by an external 
potential within a region whose size is much smaller than 
the spatial resolution of an experimental device. Then, 
the two particles would be considered as a single particle, 
which results in the overestimation of the total entropy 
production. 

The results are summarized in Fig. [21 The heat flux q 
increases from 0 to qca as K increases from zero to infin¬ 
ity. The apparent heat flux qs is different from q. When 
K is small, the hidden heat flux {qp + q-m) is negative 
so that the resolution limited measurement would over¬ 
estimate the entropy production. When K is large, the 
hidden heat flux is positive so that the resolution limited 
measurement would underestimate the entropy produc¬ 
tion. The threshold Kq is determined from the condition 
q = qs, which yields that 


Ko = 


71(71 + 72) 

2 m 


( 33 ) 


Surprisingly, the fast variable contributes to the heat flow 
by the same amount as the slow variable at all values 
of K. This hidden entropy production by the fast vari¬ 
able explains the discrepancy between qoo and grigid- The 


The conclusion on the hidden entropy production does 
not depend on the specific form of the interaction po¬ 
tential V{xi,X2)- We also considered the system with 
V{x\,X2) = ^{xi — X2)‘^ + \kx\ + \kx2- That is, each 
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particle is trapped with the additional harmonic potential 
with the coupling constant k. This case is also analyzed 
in Appendix |Al The total heat flux and the components 
are given in (IA 5 I) . With this setting, we also obtain the 
same result that q = 2grigid and that qs = qp = frigid 
with qm = 0 in the K ^ oo limit. 


V. SUMMARY AND DISCUSSIONS 

We have investigated the role of a fast variable in the 
entropy production. For the harmonic system, the total 
entropy production is separated into the contributions 
from the slow variable (center of mass coordinate) and 
the fast variable (relative coordinate), and the mixed con¬ 
tribution. We found that the hidden entropy production 
due to the fast variable remains finite even in the infinite 
coupling limit. 

Our study reveals the difference between the rigid rod 
system and the harmonic system in the infinite cou¬ 
pling limit. In the former system, the relative coordinate 
X = xi —X2 is constrained to be zero without any dynamic 
fluctuations. In the latter system, the amplitude of the 
relative motion is also negligible {x ~ How¬ 

ever, the velocity fluctuations are non-negligible, which 
results in a finite additional hidden entropy production. 

The hidden entropy production could be relevant for 
heat engines. When one models a heat engine, some de¬ 
grees of freedom may be ignored for the sake of simplicity. 
For example, in the Feynman ratchet, vanes and a ratchet 
are assumed to be linked with a rigid axle [iMl , which 
ignore relative motions of the two parts. Our results sug¬ 
gest that the entropy production in such a model should 
be less than the actual one because of the hidden entropy. 
The entropy production is the origin for the loss in the 
efficiency of a heat engine. Our result warns that all the 
fast variables should be taken into account in analyzing 
the efficiency of a heat engine. 

Appendix A: Evaluation of the integral in (I17II 

The denominator (numerator) of the integrand in ( 1171 ) 
is the 8th (6th) degree polynomial in w. Thus, one can 
add a half-circular contour of the infinite radius in the 
complex uj plane to make a closed contour as shown in 
Fig- El The contour integral is determined by the residues 
of the poles that are located at the roots of D{uj) and 
D{—u!) inside the contour. 

We first consider the roots of the quartic polynomial 
D{uj). A trivial root is located at w = 0 . The others 
are the solutions of g{uj) = m?uj^ + im(7i -|- 72)^^ — 
(7172 + 2 mK)uj — i(7i +j2)K = 0 . It is more convenient 
to work with 5(H) = ig(io = ifl) = + m(7i -I- 

72)H^ -I- (7172 + 2 mK)Q + (71 -I- 72)7^ whose coefficients 
are real and positive. Positivity of coefficients guarantees 
that real positive roots do not exist and that one of the 
roots, denoted by fli, must be real and negative. The 


Im{w} 



FIG. 3. Integration contour for the evaluation of nil in com¬ 
plex u) plane. Simple poles are marked with x symbols. The 
non-zero roots ^^4,5,6 of D{—uj) are in the upper half-plane, 
while the non-zero roots a;i,2,3 of are in the lower half¬ 

plane. 


remaining two roots, denoted by Q .2 and O 3 , may be real 
and negative, or a pair of complex conjugates. In order to 
determine whether they are inside the contour in Fig. |31 
we need to know the sign of the real parts of H 2 and 
H 3 . The three roots of g(H) satisfy Hi -|- H 2 -I- O 3 = 
-(71 + 72 )/m. Note that g(-(7i -h 12 )/m) = -(7172 + 
-I- 72 )/m < 0 and that g(x) > 0 for any real x > 
0. This implies that Hi > —(71 +^ 2 )l'ni and Re{H 2 , 3 } < 
0. Consequently, the nonzero roots Wfc = iHfc (k = 1, 2, 3) 
of are located in the lower half-plane. The nonzero 
roots oj^+k = —Wfc {k = 1,2,3) of D{—ui) are located in 
the upper half-plane (see Fig. [S]). 

We identihed all the roots of D{uj)D{—uj). Among 
them the degenerate root a; = 0 do not contribute to 
the integral in (1171) because it is canceled with in the 
numerator of the integrand. The contour in Fig. [3] in¬ 
cludes only UJ 4 , UJ 5 , and uje- Therefore, the integral in 
dni is equal to 

Mij = 2Tri [Res(a; 4 ) -|- Res(w 5 ) -I- Res(a; 6 )], (Al) 

where Res(a;fc) denotes the residue of the integrand func¬ 
tion in (El at the simple pole The roots and the 
residues are easily evaluated by using a software capable 
of symbolic algebra, which yields the results in (flSl) . 

We can generalize the analysis to the system with the 
potential energy 

V{xi,X 2 ) = ^{xl + xl) + Y{xi-x- 2 f (A2) 

instead of the potential energy in ([3|. For simplicity, we 
set 7i = 72 = 7. This system has the the modified matrix 
F(w) whose elements are Fii(a;) = ^ 22 ( 0 ;) = —— 
ijoj + k + K and Fi 2 {uj) = F 2 i{uj) = —K. Then, D{uj) = 
det F(a;) has four roots at 


a;i,2 = [—*7± V^mk — 7^]/(2m), 

W37 = [—17 ± y^ 4 m(fc -I- 2 K) — 7^]/ { 2 m) 
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in the lower half-plane, and D{—uj) has four roots at 
u! 4 +j = —ujj (j = I,-- - ,4) in the upper half-plane. 
Therefore, the integral corresponding to dig is deter¬ 
mined by the residues at The results are 


In the K ^ oo limit, we obtain that qm —> 0 and qs = 

qp ^ ^oo/2 = Qrigid- 


Mn = 


m 


M 22 —- 

m 

Mi 2 = M21 


K\T4-T2) 

' 2{mK^ + + K)y 

if2(ri -T 2 ) 

’ 2 {mK^ + + K)y 

= 0 . 


(A4) 


Inserting these into Eqs. (I12|) . (|27)) . (|28)) . and (l29l) . we ACKNOWLEDGMENTS 

obtain that 


2{mK^+ + 

qs = qp = t^(Ti - T2), 

4m 

7 y{k + K) _ 

2m {mK'^ + 'y'^{k + A^)}— - -^ 
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